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Abstract 



We introduce the phase model on a lattice and solve it using the alge- 
braic Bethe ansatz. Time-dependent temperature correlation functions of 
phase operators and the "darkness formation probability" are calculated 
in the thermodynamical limit. These results can be used to construct 
^ , integrable equations for the correlation functions and to calculate their 

asymptotics. 
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Phase operators were intensively studied in quantum optics g, They 
can be defined by the following commutation relations 

[7V,0+]=0+, [TV, 0] = -</., =7r, (1) 

where tt is the vacuum projector 

^=(|0)(0|) 

and N is the number of particles operator. 

The phase model is a model of interacting phase operators on a lattice. It 
was constructed in Q as a limit case of the q-boson hopping model j|, |^ . The 
Hamiltonian of the phase model has the following form 

1 

^=-2 E('^"'^"+i + - 2A^«), (2) 

where operators (pj , (j)^ , Nj commute in the different sites and satisfy the com- 
mutation relations (|l|) in the same site. 

The complete set of eigenvectors for the model can be obtained by means 
of the algebraic Bethe ansatz (see |Q] and references therein), i-operator of the 
model has the form 

This operator satisfies the bilinear relation 

R{p, s)Ln{p) i„(s) = L„(s) ® Ln{p)R{p, s), 

in which R(jj,q) is the 4x4 matrix i?-matrix. The non-zero elements of this 
i?-matrix are 

Rll{P,s) = R44{P,S) = f{s,p), 

R22{p,s) =i?33(p,s) ^g{s,p), 

R23{P,s) = 1 

and 

fip, s) = i I ; g{p, s) = . ' . (4) 
2sm(-2^j 2sm(^^j 

The Bethe equations for the model 

N 

exp{i{M + N)pj} = {-lf-'exp{iY,Pk} (5) 

fe=i 

(j = 1, A^) are exactly solvable: 

M + N ' ^' 
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where Ij are integers or half-integers depending on N being odd or even. 
The iV-particle eigenenergies of the Hamiltonian — H — flN (||) are 



N 



EM^Y.^h{pk)-fl); Hp) = 2sin2(p/2). (7) 



fc=i 



Here p, is the chemical potential, < /2 < 1. 

The thermodynamics of the model we shall consider for the case when the 
total momentum P = X^I^i Pk is zero P — 0. 

The thermodynamics of the model is handled in the standard way. The 
ground state energy of the model at finite temperatures /3~^ is determined 
through the solution of the nonlinear integral equations 

TT 

e{p) = h{p) - fi - (2^/3)-i / ln(l + e^<P^)dp, (8) 



2Trp{p){l + e^<P^) = 1+ / Pip)dp. 



The function p{p) is a quasi-particle density while e{p) is the excitation energy. 
The pressure is then 

V = (27r/3)-i J ln(l + e'^'^P'^)dp (9) 

— TT 

and the density is 

TT 

dV f 

D = — = J p{p)dp. (10) 

— TT 

So we have 

<p) = Hp) ~fi-v (11) 

and the quasi-particle density has the Fermi-like distribution 

27rp{p) = {l+D){l + e^''^PY\ (12) 

At zero temperature — the ground state is the Fermi sphere —A < p < 
A (A < tt) filled by the particles with the negative energies £o{p). The pressure 
and density are now 

A A 

Po = -(27r)-i J eaip)dp, Do ^ J Po{p)dp. 
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From (|Tl]) and (|T|) we have 



eo{p) = h{p)-fl~Vo, eo(±A)=0; 

2Trpo{p) - (l + Do). 

From the last equation we can express the bare Fermi momentum A as the 
function of density 



A 

The Fermi velocity v is equal: 



1 + ^0 



-d + Z^orsin.'^^" 



2^po(A) l + Do' 

When A ^ 0(/[t — > 0), Dq — > and Pq ^ as must be expected. When 
A —>■ t: (p ~* 1) all the vacancies are now occupied by particles Dq —>■ oo, Vo ^ 
and the model (||) is the classical XY chain in this limit |^ , . 

The correlation functions for the phase model can be represented as Fred- 
holm determinants of integral operators. The similar representations were ob- 
tained recently for the model of impenetrable bosons [|[ |l^ and for the XXO 
Heisenberg chain . It was shown that such results can be used for calculation 
of the asymptotics of correlation functions. 

The simplest correlation function is the emptiness formation probability. It 
can be defined as a probability of the states such that there are no particles in 
the first m sites of the lattice (see e.g. Q]). We will call it darkness formation 
probability in the case of the phase model. Using the algebraic Bethe ansatz 
and the representation for the scalar products of the Bethe states |^ one can 
represent this function as a Fredholm determinant 

T{m,P) = {1 + D)det{i - M), (13) 

where / is the identity operator and M is an integral operator 

{Mf){p) = J M{p,q)fiq)dq, 

— TT 

with the kernel 

M(p, g) ^ f V^ f" ~ f V^), (14) 

iTT sm ^(P — q) 

where v(j), P) = (1 -I- exp(/3e(p)))~^. 

It can be easily shown using the representation ( p^ ) that T(m, P) < 1. 
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At zero temperature one has 

To{m) = (1 + Do)detii - Mq), (15) 
where Mq is an integral operator 

A 

(Mo/)(p) = y Mo{p,q)f{q)dq, (16) 
-A 

with the kernel 

Mo{p,q)^- . (17) 

The time-dependent correlation function of phase operators can be also ex- 
pressed as a Fredholm determinant. This representation can be obtained using 
the representation for the form- factor p^ . We will consider the temperature 
mean values 

f+{p,m,t)^{4>„.+i{t)^tm0. (18) 
r(/3,m,t) = (0++i(t)^i(O))^, (19) 

where 

(t>m{t) = exp[ii7^t]0„exp[-ii7^,<]. 
This correlation function can be represented in the following form 

oo ^ 

/W(m,i,/3) =cxp(-^) I e'^™-')^/^^^)(M,Ae)de. (20) 

The functions /i(±'(/, /3, 8) can be written as a Fredholm determinant 

/iW(?,i,/3,e)- [G{l,t) + ^ Aei{i + V T xR±)U=o, (21) 
where V and R are integral operators 

(V^/)(P)= 2^ / V{p,q)f{q)dq, 

(22) 

(^"±/)(p) = 2^ / R^{p,q)f{q)dq, 

— 7r 

with kernels 

sm^(p — 2 
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- ^{E+il, t, q, /3, e) + E^{1, t, q, /3, Q))E^il, t,p, P)), (23) 

R+ {p, q) ^E+{1, t, p, f3, e)E^ {l,t,q,p,e), (24) 

R- (p, q) =E^{1, t, p, P)E^ (l, t, q, /?), (25) 

where the functions G{1, t),E^{l, t,p, 0),E^{1, t,p, /?, 6) and E^{1, t,p, (3, 6) are 
defined as 

G{l,t) = 2^ y exp{ilq ~ ite{q))dq, (26) 

— TT 

£;(M,P,e) = fv.p. / ^J^^^p^dq + t^.nhp-Q)eM^lp-^te{p)), 
2tt J tani(g-p) 2 

(27) 

E.{l,t,p,(3) = V^P^exp (^-iLp + iie{p)^ , (28) 
i?+(Z,t,p,/3,e) = ii?_(/,t,p,/3)x 

X ((£;(?, t,p,e) + e-'Q£;(; + i,t,p,e)) + i(G(/,t) + e-'®G(; + f,i))), ^ ' 

E-{l,t,p,l3,Q)^ ^E^il,t,p,p)x 
X ((£:(Z,t,p,e) + e'®£;(; - l,t,p,e)) - i{G{l,t) + e'^^GH - l,t))) . ^ ^ 



Although the kernel (|2^) seems to be very complicated it has the form which 
allows to use the method proposed in for calculation of the asymptotics. 
In the case t = 0,m > 1 the correlators have the following form 

/W(m,0,/3) = exp(^)i±^^ /e^(™-')«/.W(^0,/3,e)de, (31) 

I— in „ 

— TT 

h^^\l, 0, 13, e) = ^det(/ - {) + xf±)|,=o, (32) 
where the integral operators v and possess the kernels 

-sin ^(p -q)+ cxp[i(e - sin ^(p - q) 



v{p,q) = ^Mp) '-^^—^ ^ (33) 

sm2(p-'7) 

r"(p,'?) = V^e-*^(f+«)VK9)- (35) 

Comparing the results obtained in this paper with the representations for 
the correlation functions of the other models (||^, |l^, 0) one can see that these 
representations are rather different and there are new difficulties in our formulae. 
The models considered before were the free fermion points of the more general 
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integrable models. The phase model being a limiting case of the g-boson hopping 
model is not a free fermion point of this model and this can explain the new 
difficulties appeared in the correlation functions. The very peculiar property 
of the correlation functions of the phase model is the dependence on the full 
density of particles. 

This work was supported by the grant number 95-01-00476a of the Russian 
Foundation of Fundamental Research. A.I. and N.K. thank Ecole Normale Su- 
perieure de Lyon for hospitality. N.B. would like to acknowledge the hospitality 
at the Research Institute for Theoretical Physics University of Helsinki. 
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